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ABSTRACT 

We analyse the polarized emission at 1.4 GHz in a 3° x 3° area at high Galactic 
latitude {b ~ —40°). The region, centred in (a = 5^, (5 = —49°), was observed with 
the Australia Telescope Compact Array radio- interferometer, whose 3-30 arcmin an- 
gular sensitivity range allows the study of scales appropriate for Cosmic Microwave 
Background Polarization (CMBP) investigations. The angular behavior of the diffuse 
emission is analysed through the E- and B-mode angular power spectra. These follow 
a power law Cf ex with slopes = -1.97 ± 0.08 and /^^ = -1.98 ± 0.07. The 
emission is found to be about a factor 25 fainter than in Galactic plane regions. The 
comparison of the power spectra with other surveys indicates that this area is inter- 
mediate between strong and negligible Faraday rotation effects. A similar conclusion 
can be reached by analysing both the frequency and Galactic latitude behaviors of 
the diffuse Galactic emission of the 408-1411 MHz Leiden survey data. We present 
an analysis of the Faraday rotation effects on the polarized power spectra, and find 
that the observed power spectra can be enhanced by a transfer of power from large 
to small angular scales. The extrapolation of the spectra to 32 and 90 GHz of the 
CMB window suggests that Galactic synchrotron emission leaves the CMBP £^-mode 
uncontaminated at 32 GHz. The level of the contamination at 90 GHz is expected 
to be more than 4 orders of magnitude below the CMBP spectrum. Extrapolating to 
the relevant angular scales, this region also appears adequate for investigation of the 
CMBP 5-modes for models with tensor-to-scalar fluctuation power ratio T/S > 0.01. 
We also identify polarized point sources in the field, providing a 9 object list which is 
complete down to the polarized flux limit of Sf^^ = 2 mJy. 

Key words: cosmology: cosmic microwave background - polarization - radio con- 
tinuum: ISM - diffuse radiation - radiation mechanisms: non-thermal. 



1 INTRODUCTION 

Polarized Galactic synchrotron emission is one of the most 
important foregrounds in measuring the Cosmic Microwave 
Background Polarization (CMBP) up to about 100 GHz. 
Above that frequency dust emission becomes the dominant 
contaminant. Moreover, as the fractional polarization of the 
synchrotron emission can be 30% and as CMBP is just few 
percent of the CMB anisotropy, the degree of contamina- 
tion is expected to be worse than in temperature anisotropy 
measurements. 

The study of this foreground is thus crucial for CMBP 
experiments. It will allow an estimate of the contamina- 
tion level, and will aid developing cleaning procedures to 
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remove its contribu tion from the cosmic signal (e.g. see 
iTegmark et aLlEoOOl and references therein) . 

The CMBP emission peaks on sub-degree 
ang:u lar scales in th e 5-30 arcmin range (e.g. 
see IZaldarriaga. Soergel SeliakI IIQQTT) . Consequently, 
observations of CMBP can be carried on in small sky 
patches that are large enough to allow good statistics 
(5°-10° wide) but, at the same time, small enough to 
minimize environmental systematics, e.g. spillover from 
ground emission. The latter could swamp the detection of 
the faint (few /xK) cosmic signal. Small regions, moreover, 
have the additional advantage of allowing a part of the sky 
low in synchrotron contamination to be selected. 

The synchrotron emission can be best observed at low 
frequency, where it dominates other diffuse components 
(dust, free-free, and CMBP itself) and where the signal 
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strength allows easier detection by radiotelescopes. To date, 
analysis of the Galactic synchrotron has been performed 
by several surveys at frequencies up to 2.7 GHz. These 
have been mainly concentrated on the Galactic plane, where 
the sign al is greater. The So uthern Galactic Plane Survey 
(SGPS, iGaensler et al.l l2Q0lh and the Canadian Galactic 
Plane Survev fCGPS. lTavlor et al .1 2003 ^ are interferometric 
surveys whose main goal is to map almost all the Galactic 
plane at 1.4 GHz down to a 1 arcmin resolution. These sur- 
veys provide deep insight into many effects typical of Galac- 
tic polarized emission (e.g., polarization horizons, Faraday 
screens). However, because of their interferometric origin, 
these surveys are sensitive to angular scales no larger than 
the 30 arcmin of the telescope primary beam. 

The largest angular scales are instead covered by 
a single dish project: the Effelsbe rg 1.4 GHz M e dium 
Galactic Latitude Survey (EMLS, lUvaniker et al.1 Il999l 
and iReich et al.l l2004l) . This covers the Galactic plane up 
to medium latitudes < 20°). The approximately 10 ar- 
cmin resolution allows EMLS to overlap the scales accessible 
by the two previous interferometric surveys, and their com- 
bination will provide full information on the Galactic plane 
down to 1 arcmin. 

The G alactic plane h as been also surveyed at higher fre - 
quency bv lPuncan et al.l (1997) and iDuncan et al.1 (119991) . 
They covered about a half of the plane at 2.4 and 2.7 GHz, 
respectively, with the latitude coverage extends to \b\ < 5°. 

The mid Galactic latitudes has been par- 
tiall y surveyed at 350 MHz by the W esterbork sur- 
vey (jHaverkorn. Katgert de Bruvnll2003r) . which mapped 
the polarized emission at longitudes / ^ 140°-170° and up 
to 5 = 30°. Because of their lower frequency, these data are 
likely to be more affected by Faraday rotation effects. 

The analysis of the Galactic plane region has provided 
the first information on the angular behavi or of t he Galac- 
tic synchrotron radiation fsee lBruscoli et aI1l2002l and refer- 
ences therein). However, the optimal locations for observa- 
tions of the CMBP are at high latitudes, where the emission 
is low. High latitude regions are not as well studied, although 
more recent surveys are starting to fill this gap, at least at 
1.4 GHz. 

Sparse observations at all Galactic latitudes were 
made in 1960s and 1970 s (e.g. iBaker W ilkinson, Il974l: 
iBrouw Spoelstral Il976l . and referenc es therein). Among 
them, the observations presented by iBrouw Spoelstral 
(|l976l) covered all Galactic latitudes at 5 frequencies be- 
tween 408 and 1411 MHz. For decades these have repre- 
sented the largest data set of polarized measurements out of 
the Galactic plane. Although undersampled and not suitable 
for a full investigation of the synchrotron characteristics (es- 
pecially when dealing with the angular behavior) some anal- 
yses of the angular power spectra have been carried out in 
the best sampled portions (Bruscoli et al. 2002). Moreover, 
as we will see in this work, this data-set can provide useful 
information on the frequency behavior of the large angular 
scale structure. 

The whole of the North ern and Southern sky 
are being mapped by Wolleben et al.l (120041) and 
iTestori. Reich ReichI (l2004l) . respectivelv. at 1.4 GHz. 
These surveys are very precious, since they will provide the 
first all-sky map of the polarized emission at 1.4 GHz down 
to the degree scale. They are sensitive to the large angular 



scale structure of the polarized emission in interesting (for 
CMBP purposes) low emission regions. However, their 
angular resolution (larger than 30 arcmin) and sensitivity 
(about 15 mK) do not allow the analysis of the faint areas 
required by sub-degree CMBP experiments, which require 
better resolution (finer than 30 arcmin) a nd sensitivity 
(rms signal about 10 mK, see lBernardi et al ]|2003). 

All of the aforementioned surveys are at low frequencies, 
where the Faraday rotation effects can still play a significant 
role. In fact, Faraday rotation can introduce a randomization 
of polarization angles which transfers power from large to 
small angular scales by modifying the polarized emission 
pattern. Although this transfer of power has been claimed 
to explain some results at frequencies up to 1.4 GHz (e.g 
iTucci et alll2002^ . it has not been quantitatively studied. A 
detailed study is needed to evaluate its effect before we can 
safely extrapolate up to the mm-wave cosmological window. 

The patch of sky near a — 5^ and 6 — —50° has been 
identified as an interesting area for CMBP investigations: it 
is at high Galactic latitude {b ^ —40°) and, from an analy- 
sis of the Rhodes/HartRAO 2326-MHz radio continuum sur- 
vey ( Jonas, Baart & Nicolson 1998), it is expected to have 
a synchro tron emission at the low level requi red for CMBP 
studies (.BernardilEooi : ICarretti et al.ll2002^ . In fact, this 
region was chosen for the total intensity observation by the 
BOOMERanG experim ent and h as been selected as target 
for Ba R-SPOrt (Cortiglioni et al.k2003 ) and BOOMERanG- 
B2K dMasi et al ll20o3v ! 

This patch has been observed in polarization at 1.4 GHz 
to a sensitivity t o allow detection o f the polarized syn- 
chrotron emission ( Bernardi et al.l2003 ). This represents the 
first detection of this emission at high Galactic latitude at 
this frequency and in the angular scale range useful for 
CMBP analyses (5-30 arcmin). 

These observations, along with an initial analysis of 
the implication s for C MBP studies, have been presented by 
iBernardi et al.l (|2003l) . In the present paper we perform a 
more detailed analysis of the polarized emission in this area 
and reach firmer conclusions regarding CMBP implications. 

We quantatively analyse the effects of the polarization 
angle randomization introduced by Faraday screens, find- 
ing that a power transfer from large to small angular scales 
occurs. This results both in a steepening of the (angular) 
spectral index and an enhancement of the observed emis- 
sion on the angular scales our observations are sensitive to 
(i.e. 3-30 arcmin). This implies that our measurements rep- 
resent an upper limit of the intrinsic polarized emission. 

We present the angular power spectra of the E- and 
B-modes. Comparison of these with other surveys suggests 
that the observed patch is an intermediate state between 
significant and negligible Faraday rotation effects. 

Finally, we extrapolate the spectra up to the frequency 
range of CMB measurements (30-90 GHz). This suggests 
that the Galactic synchrotron emission is low enough in this 
area that it should not be an issue for both the £^-mode 
(from 30 GHz) and the fainter B-mode (provided a tensor- 
to-scalar power ratio T/S > 0.01). 

The paper is organized as follows. In Section |21 we 
present a quick summary of the observations. In Sections El 
and ^ we analyse the point sources detected in the field and 
the Rotation Measure (RM). Section [^discusses the power 
transfer from large to small angular scales as a result of 
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randomization of polarization angles. It considers how this 
changes the polarized power spectrum. The E- and B-mode 
power spectra of the synchrotron emission in the patch will 
be presented in Section |S| together with an analysis of the 
role of Faraday rotation effects. The implications for CMBP 
investigations are discussed in Section Q Summary and con- 
clusions are finally presented in Section |H| 



2 OBSERVATIONS 

The region centred on {a - 5^, ^ — has 
been observed at 1.4 GHz by iBernardi et al.l (|2003f) 
using t he Australia Telescope Compact Array (ATCA, 
iFrater. Brooks & Whiteoak 1992|), an East- West synthesis 
interferometer situated near Narrabri (NSW, Australia), op- 
erated by CSIRO-ATNF. 

The observation was performed as a 49 pointing mo- 
saic covering a 3° x 3° region. The array configuration used 
was the so-called EW214 including spacing from 30 m to 
214 m. This provides sensitivity on angular scales ranging 
from ^ 30 arcmin down to the resolution of ^ 3.4 arcmin. 
The system provides the four Stokes parameters /, Q, U, 
V. Details of the observations are listed in Tableland pre- 
sented in Bernardi et al. (2003). 

The maps of the Stokes paramete rs /, Q a nd U, along 
with the polarized intensit y — a/Q^ + IP , a re shown 
in Figure Differing from IBernardi et all (|2003l) . here no 
subtraction of polarized point sources has been performed. 



Table 1. Main characteristics of the 1.4 GHz observations. 



Central Frequency 


1380 MHz 


Effective Bandwidth 


205 MHz 


Array Configuration 


EW214 


Angular Sensitivity Range 


3.4-30 arcmin 


Location (J2000) 


a = h^,5 = -49° 


Area Size 


3° X 3° 


Observation Period 


June 2002 


Effective Observing Time 


~ 70 h 


Sensitivity (flux) 


0.18 mJy beam"-*^ 


Sensitivity (temperature) 


3.2 mK beam-i 



The polarized emission P is patchy in nature, and dis- 
tributed over the whole field. An exception is a bright fea- 
ture in the South- West corner of the area: it is more fila- 
mental, being about 1° long. The Stokes / image does not 
show any particular diffuse structure at the same coordi- 
nates, although the point source contamination makes the 
comparison hard. As in other high resolution observations 
fe.g. IWieringa et al.lll993l : iGaensler et al.ll20(Hl) . a Faraday 
screen is likely acting along the line of sight. This acts as a 
small scale modulation of a relatively uniform background, 
generating the apparent filamental structure on small angu- 
lar scales. 

The mean pol arized emission Prms = ^y{Q^) + (lP) — 
11.6 mK found bv IBernardi et all (120031 ) is well above the 
beam sensitivity (S/N ~ 3.5) allowing an analysis of the 
synchrotron emission features in the area. 
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Table 2. Position, total intensity (/) and polarized intensity (/^) of the 9 sources detected in the field. The rms-error is the same for 
both the two intensities and corresponds to the beam-sensitivity (0.18 mJy beam"-*^). Polarization angle (^) and polarization degree 
(n = /I) are also reported. Counterparts were located in radio catalogues. Where no radio-counterpart was found, the closest optical 
object found in both the APM and 2-MASX surveys is listed (see text for details). The last column provides the distance from the 
counterpart. 



source 


RA 
J2000 




DEC 
J2000 




/ 

[mJy] 


IP 
[mJy] 






n 


counterpart 


distance 
[arcmin] 


1 


5^01^01^ 


.4 


-48°31'03'' 


.0 


164.10 


10.72 


15.4° 


±0.5° 


6.5% 


PMN J0501-4831 


0.5 


2 


5^05™58^ 


.8 


-49° 11^43'' 


.0 


98.52 


6.63 


21.5° 


±0.8° 


6.7% 


2MASX J05055521-4910485 


1.1 


3 


4^58^34^ 


.5 


-48° 34^40'' 


.2 


130.20 


3.71 


43.8° 


± 1.4° 


2.9% 


APMUKS(BJ) B045704.24-483820.7 


1.8 


4 


5^00™38^ 


.2 


-49° 12^29'' 


.5 


139.70 


3.40 


43.6° 


± 1.5° 


2.4% 


PMN J0500-4912 


0.2 


5 


4^59^37^ 


.8 


-48°43'13'' 


.6 


38.21 


2.78 


41.3° 


±1.9° 


7.3% 


APMUKS(BJ) B045822. 87-484738 


0.9 


6 


5^07^08^ 


.8 


-49°29'32'' 


.0 


20.73 


2.75 


-27.5^ 


' ±1.9° 


13.3% 


2MASX J05070829-4929268 


0.1 


7 


4^57^17^ 


.7 


-48° 34^58'' 


.5 


23.93 


2.32 


42° 


±3° 


9.7% 


APMUKS(BJ) B045558. 02-483812. 5 


1.3 


8 


5^07™24^ 


.5 


-50° 13^14'' 


.4 


123.70 


2.18 


40° 


±2° 


1.8% 


PKS 0506-502 


0.2 


9 


4^59^29^ 


.7 


-48°31'56'' 


.9 


55.48 


2.01 


44° 


±3° 


3.6% 


PMN J0500-4912 


0.3 



3 POINT SOURCES 

Even though the field has been selected to minimize their 
contamination, a few point sources are present in the po- 
larized intensity map (Figure Q). To identify them, we 
fit the maxima/minima in Q and U maps with a 2D- 
Gaussian beam and a constant value using MIRIAD's task 
IMFIT (iSauTt. Teuben S^, Wrightlll99,^ . The constant value 
is adopted to account for the background emission. We 
keep only those sources where the resultant Full Width at 
Half Maximum (FWHM) is compatible with the synthe- 
sized beam. Sources providing larger FWHM are instead 
discarded, the fit being not able to separate the background 
emission. 

For the detections, we retain only those point sources 
where the emission is 3 times the rms signal of the dif- 
fuse component {S/Prms > 3). This gives good confidence 
that we avoid confusion between point sources and statis- 
tical fluctuations of the diffuse background, while ensuring 
a solid detection above the instrumental noise {S/N > 10). 
Moreover, this criterium allows us to provide a complete list 
down to the polarized flux of S^^^ = 2.0 mJy. 

Given to these criteria, we find 9 sources: these are listed 
in Table 

We have attempted to identify the sources in pub- 
lished catalogues such as the Parkes-MIT-NRAO (PMN, 
iGriffith X. Wr ightl ll99,'j^ survey. PMN was performed at 
4.85 GHz, so that, assuming a mean spectral index a = 
0.7, its S4.85 = 42 mJy flux- limit corresponds to about 
Si A 100 mJy at 1.4 GHz. Accepting this, we would ex- 
pect sources 1, 2, 3, 4 and 8 to be present in the PMN 
catalogue. However only sources 1, 4 and 8 were found: we 
are not able to find any identification for sources 2 and 3. 
It is worth noting that their fluxes are near the catalogue 
limit: a steeper spectrum would result in PMN failing to 
detect them. Conversely, source 9 has a PMN counterpart, 
even if the measured 1.4 GHz flux is about a half of the limit 
extrapolated from 4.85 GHz. 

When a PMN counterpart could not be found, 
we looked for the clo sest objec t present in the opti- 
cal/IR catalogues A PM (|Maddox e t al. 1990) and 2-MASX 
(jHuchra et alJbOO.^) . provided it lies within a FWHM of 



the detected radio source. Possible counterparts have been 
found for all the missing objects. These are listed in Tabled 
It is worth noting that all of these sources are classified as 
normal galaxies. 

Given the limited number of sources, no general prop- 
erties can be extracted. We simply note that only one source 
exceeds the 10% polarization. 



4 ROTATION MEASURE 

The magnetic field parallel to the line of sight {B\\) changes 
the polarization angle by Faraday rotation. The variation 
A(f) — (f) — (f)Q with respect to the intrinsic polarization angle 
00 is given by the formula 

A0 = iWA^ (1) 

where A is the wavelength of the radiation and RM is the 
Rotation Measure. The latter is defined by the integral along 
the line of sight 

iW = 0.81 J B\\{1) ne(l) dl (2) 

which depends on B\\{1) and the free-electron density ne{l) 
in the Interstellar Medium (ISM) at various distances / from 
the observer. 

Beside changing in the polarization vector direction 
and, in turn, in the naive estimate of the magnetic field ori- 
entation, Faraday rotation can induce both depolarization 
effects and Faraday screen modulations. These can signifi- 
cantly modify the angular power distribution and pattern 
of the polarized emission. EM estimates in the patch of sky 
being consider are thus important in understanding the sig- 
nificance of these effects. 

Because the ATCA correlator produces a number of fre- 
quency channels across the observed bandwidth, we are able 
to determine values of RM from our data set. Faraday rota- 
tion has been evaluated by grouping the 26 useful channels 
of our observations in four sub-bands to form four maps of 
Q, U and V at the different frequencies of 1316, 1368, 1404 
and 1454 MHz. For each frequency, the Stokes V map sets 
the noise level for the corresponding Q and U maps. Using 
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filament structure. This is thus characterized by EM values 
larger than typical of this patch. This might be expected if 
the filament is caused by a Faraday screen. The other mea- 
surements of EM lie near the first peak. Thus 20 rad m~^ 
can be considered a typical value of the EM of the successful 
pixels. 

This result is in good agreement with other estimates 
of EM in this region. The all-sky cata log!:ue of rotation mea- 
sures for extragalactic sources by Br oten. MacLeod Valle3 
(1988) contains four sources within 5° from the centre 
of the area. These have EM values of 13, 28, 34 and 
53 ra d m~^. The rotation measure maps oflHan et al.l (|l997f) 
and J ohnston-Hollitt. Hollitt EkerJ (|2004ri are indicative 
of large-scale EM values. In our ar ea, they su gg!:est values of 
approximately 20-30 rad m" rsee lHanl[2003 for a review). 

Given these various estimates of EM, we can consider 
EM — 50 rad m~^ as a reasonable upper limit for the typi- 



(rad/m^2) 



cal rot ation measure within the patch. As in lBernardi et al.l 
(|2003f) . this value in combination with the total bandwidth 
of the ATCA observations implies a depolarization factor of 
D ^ 92%. This has a marginal impact on our main goal of 
estimating the level of polarized synchrotron emission in the 
area. 



Figure 2. Map of RMs measured in the patch. 
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Figure 3. Distributions of the measured RMs. 

MIRIAD's task IMPOL, we have obtained maps of polarized 
intensity and polarization angle at each frequency. EM val- 
ues are then determined using MIRIAD's task IMRM, that 
fits a EM map from polarization angles at the different fre- 
quencies. This task also tries to solve the Ntv ambiguity in 
the EM determination. Pixels with an error in polarization 
angle greater than 20° are discarded. 

The EM values that were successful determined by this 
process are shown in Figure |2l The coverage of the mea- 
surements in the patch is sparse, and so they cannot be 
considered to be fully representative of the area. However, 
they are scattered over all the field and thus give a useful 
indication of the broader trend. The distribution, shown in 
Figure 13 presents two peaks: the first near 20 rad m~^ and 
a second near 80 rad m~^. Apart from a small number of 
points, the measurements of the second peak relate to the 



5 POWER TRANSFER FROM LARGE TO 
SMALL SCALES 

As n oted by several aut hors (e.g. see IWieringa et a'llll993l 
and iGaensler etHIEoO lV Faraday screens can introduce a 
small scale modulation of a relatively uniform background. 
If a complex Faraday screen is interposed between the ob- 
server and a uniform polarized emission, the EM pattern of 
the screen will cause modulation of polarization angle, re- 
sulting in variable patterns of Q and U. This causes transfer- 
ring power from large to small angular scales in polarization 
maps, and generates false small scale structures. 

Figure 31 demonstrates this point using a simple model: 
the intrinsic emission is uniform, while the observed one is 
sinusoidal. This is caused by a EM which varies linearly 
with angle. A result of the process is that angular power 
that had corresponded to the large scale structure has been 
transferred down to the modulation scale. 

Overall the effect randomizes the polarized emission 
larger than a particular angular scale. 

Although this effect has been invoked as a qualitative 
explanation of structures observed, we are not aware of a 
quantitative analysis having been performed. 

Here we simulate some realistic examples to estimate 
importance and main features of the effect. First, we gen- 
erate Q and U maps starting from E- and B-mode angular 
power spectra which we assume have power law behavior 

(X £^^, 

(X (3) 

with the spectral indeces which can be 

/3x = 0.0, -1.5, -3.0, X = E,B. (4) 

These have been selected in a range wide enough to cover 
all of the cases observed to date for the polarized Galactic 
synchrotron emission and are centred on the most common 
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Figure 4. A simple model showing how a variable RM pattern 
can transform uniform emission to one modulated on smaller an- 
gular scales. The top panel shows a background with uniform 
emission on a 10° box-size (emission of Stokes Q only is assumed, 
i.e. 00 = 0°). The mid panel shows RM varying linearly. The ob- 
served Q is shown in the bottom panel, assuming a 300 MHz fre- 
quency. The observed emission is modulated on 1° scale, smaller 
than the intrinsic emission. 



Figure 5. CIC scheme for regularly gridded data. The four grid 
points ii, 22, 13 and 24 are the nearest to the pixel j where the 
quantities sampled on the grid points have to be interpolated. 
Aif^ with k = 1,4 are the opposite sub-cells used in the weighted 
mean. 



measured values ((3x 
iGiardino erallEo02l . 

We use the synf ast procedu re of the HEALPix package 
(|G6rski. Hivon Wandeltl 19991) to generate the maps given 
the power spectra. We adopt an angular resolution of about 
2 arcmin (HEALPix's parameter A^side = 2048) to match the 
angular resolution of the observed maps (about 3.4 arcmin). 
Finally, 20° x 20° square maps are extracted. 

The area we observed shows a polarization an- 
gle pattern uniform at least up to 10-15 arcmin scale 
(see I 



iBernardi et al.l 1200311 . In this simulation, we assume 
that the polarization angles are totally random above 15 
arcmin angular scale size. This is done by placing a 15 ar- 
cmin X 15 arcmin grid on the simulated map, and then 



assigning a random polarization rotation angle to each grid 
point. 

To associate a random angle with each pixel of the map, 
the random angles of the nearest grid points are linearly in- 
terpolated to the pixel position by using a Cloud-in-Cell 
(CIC) scheme (iHocknev Eastwoodl[l981ll . In order to ex- 
plain the procedure, let us define Oi as the random angle 
associated to the z^^ point of the grid and j the j^^ pixel of 
the map where the angles have to be interpolated. Referring 
to Figure [51 the CIC method consists first in finding the 4 
grid points ii-i^ which are the nearest to the pixel j and 
which define the Cell where j itself is located. The linear 
interpolation is performed assigning to j a weighted mean 
of the values sampled on these 4 points: the pixel j divides 
the Cell in four sub-cells and the weight of a grid point is 
proportional to the area of the sub-cell opposite to the point 
itself: 



y4 



.1 ^ik 



(5) 



where Oj is the quantity to be estimated at the pixel j, ik 
with /c = 1,4 are the four nearest grid points and Aij^ are 
the areas of the corresponding sub-cells. 

We have adopted this linear interpolation scheme to 
avoid the discontinuities that a simple Nearest Grid Point 
(NGP) scheme would produce in the polarization angle map. 

Figure El shows the power spectra before and after 
the randomization procedure. The results for are similar. 
For the case (Se = 0.0, there is no change in the power at any 
scale. This is expected as this case corresponds to a totally 
random pattern of Q and U (Ci = constant is the case of 
pure white noise): adding further randomization does not 
change the statistics of the data. 

In the other two cases the randomization reduces the 
power on the largest scales, while increasing it on the small- 
est ones. In particular, on the scales that the interferometer 
is sensitive to (3-30 arcmin), the power is greater or equal 
to that of the input maps. 

Randomization of the large-scale polarization angle 
thus does transfer power from large to small angular scales. 
On the small scales, the angular power detected is an en- 
hancement. Consequently what we measure is an upper limit 
to the intrinsic fiuctuation. 

The angular scale which marks the transition from re- 
duction to enhancement depends on the power law index. 
However, at least for the range explored here, it is larger 
than the scale size of the randomization. Because our map 
has a randomization scale of at least 15 arcmin, the power 
detected in the 3-30 arcmin range is enhanced or at least 
not reduced. Consequently, in terms of power spectra, our 
measurements represents an upper limit of the real emission. 

From Figure El we can consider the slope of the power 
spectrum: this power spectrum is steeper in those -^-ranges 
where the power is enhanced. Indeed, apart from a transi- 
tion range centred at the randomization scale, the modified 
spectra follows a power law but with a steeper slop e. This 
is in agreement with the results of lTucci et al.l ()2002f) . They 
found spectra at 1.4 GHz that were steeper than at 2.4 GHz, 
where Faraday effects are weaker. 
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Figure 6. power spectra for the three models described in 
the text (Pe = O-O? -1-5, -3.0) before (sohd hnes) and after 
(crosses) the apphcation of the polarization angle randomization 
procedure. 



6 POWER SPECTRUM ANALYSIS 

E- and B-modes, combinations of the tensorial 2-spin quan- 
tities Q j completely describe the polarized emission 
and have the useful charac t eristic of being; sc alar (e.g. see 
IZa,lHa,rria,ga, SeliaU 11991 Izaldarria.gal 1199^. Thus, they 
allow us to describe how the polarized signal is distributed 
across angular scales by using simple scalar spherical har- 
monics. 



E mode 




1000 



Figure 7. E- (top) and B-mode (bottom) angular power spectra 
of the polarized emission at 1.4 GHz in the observed patch. 

We prefer these descriptors instead of the scalar spectra 
of Q and U because the latter depend on the orientation 
of the reference frame. Additionally, E- and B-spectra are 
the quantities predicted by cosmological models. Their use 
for Galactic work allows a direct comparison between the 
Galactic and cosmological signals and the evaluation of the 
contamination of the latter by the former. 

We have computed the E- and B-mo de spectra by us- 
ing the Fourier technique of ISeliakI (|l997ll . The results are 
shown in Figure |71 Both power spectra are well approxi- 
mated by a power law on scales smaller than about 15 ar- 
cmin (multipoles i > 800). On larger scales, there is a turn- 
over and power becomes negligible for i < 300-400. This 
corresponds to about the 30 arcmin antenna primary beam 
size. This turn-over is a result of the spatial filtering of an 
interferometer. The ATCA has little sensitivity on scales 
larger than about the FWHM of the primary beam, has 
somewhat improved sensitivity between about FWHM and 
FWHM/2, and shows full sensitivity on scales smaller than 
about FWHM/2. 

The power law behavior covers the 3.4-15 arcmin range 
of full sensitivity of the instrument. Table Elgives the results 
when we fit the functional form 

Cf=C-oo(^) , X = E,B, (6) 

in the ^-range of 800-2800. 
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Table 3. Best fit parameters of the angular power spectra of our 
1.4 GHz data. 



Spectrum 






% 


8.0 ± 0.2 


-1.97 ± 0.08 




8.0 ± 0.2 


-1.98 ± 0.07 



Table 4. Slopes (3x of angular power spectra of the Galactic syn- 
chrotron polarized emission. The table also reports the sky area 
the slope has been computed for and the observation frequency. 
SGP S-2.4 refers to the slope computed from the 2.4 GHz sur- 
vey JPuncan et al.llT'997^ in the area covered by the SGPS Test 
Region. The last column reports the references where the slopes 
have been measured. 



Area 




frequency 


Ref. 


This area 


^ -2.0 


1.4 GHz 


this work 


SGPS 


~ -2.8 


1.4 GHz 


Tucci et al. f2002) 


SGPS-2.4 


~ -1.7 


2.4 GHz 


Tucci et al. r2002) 


Galactic plane 


- -1.6 


2.4-2.7 GHz 


Bruscoli et al. f2002) 



It is interesting to compare the amplitude and spectral 
index with information from other Galactic surveys. First, 
let us consider the SGPS Test Region ( Oaensler et al.ll2nnih . 
This covers a similar area (about 4° x 6°). Like the current 
study, the SGPS Test Region survey is an interferometric 
observation at arcminute resolution and carried out at the 
same 1.4 GHz frequency. However, its region of interest is 
located on the Galactic plane. 

The first significant di fference is in the amplitude: SGPS 
has a C^oo 5 x 10"^ (|Tucci et al.ir2 nn2V corresponding 
to a signal about 25 times stronger than in the patch of the 
current study. This is expected, as the observed patch is at 
high Galactic latitudes. Furt h ermor e this is consistent with 
the result of iBernardi et al.l (|2003f) . which found that the 
mean signal of this patch is about 10 times smaller than the 
background emission in areas close to the Galactic plane of 
the 1.4 GHz EMLS survey (Uvaniker et al. 1999). However 
it is worth noting that the difference is not as large as the 
case with total intensity: in total intensity the differences be- 
tween Galactic plane and high Galactic latitudes emissions 
can be much larg er - for instance, see the full-sky surveys at 




Another area of significant difference is the slopes of 
the power spectra (Table 3|. The power laws in the ob- 
served patch have Px —2.0, whi ch is conside rably flatter 
than 2.8 for the SGPS (iTucci et al.lE002l How- 
ever it is close to the - -1.7 found at 2.4 GHz 
in t he region covered by the SGPS, bu t using the data of 
l^ncan et al. (1997) E^^^D Ilooi) . The latter agrees 
with the mean value (3^^~^'^ ^ —1.6 ± 0.2 measured at 
2.4-2.7 GHz over a large portion of the Galactic Plane 
dBruscoU et al.ll2nn2h . 

The large difference between /jt^^^" ^-^, 0x^~^'^ and 
(3%^^^ was first addressed by Tu cci et all ^2002), who ex- 
plained the steepening of the SGPS through a Faraday- 
rotation induced transfer of power from large to small scales. 
Section [51 provides a quantitative approach to this effect: 
Faraday screens can introduce randomization in the polar- 
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Figure 8. As for Figure [S] but using the modified GIG scheme 
dealing with irregularly sampled data (see text for details). 



ization angle, which actually transfer power from large to 
small angular scales and, thus, increases the power and the 
slope on scales smaller than the randomization scale size. 

Faraday rotation has a square dependence on the wave- 
length and the effects at 2.4 GHz are expected to be less 
significant than at 1.4 GHz. The flatter slope at 2.4 GHz 
can be interpreted as a reduction of the effects of Faraday 
rotation. The slopes most probably are near the intrinsic 
value. 

Moreover, in spite of a large variation of RMs across the 
Galactic plane, the slopes measured at 2.4-2.7 GHz have a 
small dispersion (A/3j^~^"^ = 0.2), indicating a good stabil- 
ity around the mean value. This further supports the idea 
that the 2.4-2.7 GHz slopes are slightly modified by Faraday 
effects and are near the intrinsic values. 

Our results are lie between the two previous values. 
In fact, pulsar and extragalactic source measurements show 
that the RM becomes less significant at high Galactic lati- 
tudes: starting from the typi cal value \EM\^ 200 rad m~^ 
at the Galactic Plane (e.g. IGaensler et al. || 200lh. through 
20-30 rad m~^ at mid Galactic latitudes (iHan et al.l I1997I: 
IJohnston-Hollitt et £111 ^004 and down t o 10-20 rad m~^ i n 
the North Galactic Pole (NGP) (e.g. see Sun Hanl l2004ri . 
Hence we would expect to see less impact of Faraday rota- 
tion in the observed patch than in the SGPS. We conclude 
that our observed patch is in an intermediate domain be- 
tween significant and negligible Faraday rotation effec ts. 

W e have analysed the data of iBrouw Spoelstral 
(Il976l), which has led us to a similar conclusion. In more de- 
tail, the data of lBrouw Spoelstral (|l976f) give polarization 
observations of about a half of the sky at five frequencies: 
408, 465, 610, 820, 1411 MHz. Having the desirable charac- 
teristics of covering all the Northern Galactic latitudes and 
a wide range in frequency, these data allow us to explore 
behavior in both frequency and latitude. An issue with this 
data-set is that it is sparsely sampled and arranged on ir- 
regular grids. While the lowest frequency case is tolerable 
(the best sampled areas have about 2° sampling versus a 
FWHM = 2.3°) the situation is worst at 1411 MHz, where 
FWHM = 0.6°. Nevertheless, regridding and smoothing on 
angular scales larger than the sampling distance allow the 
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analysis of the characteristic of the emission at least for the 
large scale structure. 

To account for irregular sampling, we regrid the data 
on regular HEALPix maps via linear interpolation, using a 
variant of the CIC scheme. The standard CIC deals with 
data sampled on a regular grid, as described in Section |5] 
On the other hand, here the data are sampled on irregular 
grids. Let (Q^, Ui) be the i^^ sample of the grid and j the j^^ 
pixel of the map to be interpolated. Referring to Figure |H| 
we generalize the CIC scheme finding the 3 grid points (ii, 
^2, is) which define the smallest triangle containing j. Con- 
necting j to ii, i2 and zs, three triangular sub-cells can be 
identified. Similarly to Eq. 0, the linear interpolation is 
performed assigning to j an average of the values sampled 
on these 3 grid points weighted for the area of the opposite 
sub- cell: 



X 



(7) 



where Xj is the quantity to be estimated at the pixel j. 

The deformations of the 2-sphere space can generate 
depolarization due to the projection of Q and U onto the 
local reference frame of parallels and meridians if a simple 
mean of the data is performed. Clos e to the NGP, t his will be 
a particular issue. Following Brusc oli et al.l (|2002l) . to avoid 
these projection effects, we perform a parallel transport of 
the polarization vectors (Qi^,, Ui,^) onto the interpolation 
point j before averaging the data. 

To account for the irregular sampling distance, the 
interpolated data are smoothed with a Gaussian filter of 
4° FWHM, leading to maps able to describe the large scale 
distribution of the polarized emission. 

The resultant maps at the 5 frequencies are shown in 
Figure El A feature is clearly visible in the Fan region at 
Galactic longitude I ~ 150° at all frequencies. This region 
is close to the area where the line of sight is nearly per- 
pendicular to the local Galactic magnetic field. Hence, there 
is only a small parallel magnetic field component, and so 
small Faraday rotation effects. Consequently the presence 
of a large non- depolarized region is not surprising, even at 
the lowest frequency. 

For the current study, the area of most interest is at 
very high Galactic latitudes. The North Galactic Spur, while 
depolarized at 408 MHz, becomes evident with increasing 
frequency. Moreover, this structure is apparent at the high- 
est latitudes first and propagates toward lower ones as fre- 
quency increases: while at 610 MHz the large scale polarized 
structure is evident only close to the NGP, at 820 MHz it is 
present down to 6 ~ 60° and reaches 6^40° at 1411 MHz. 

The polarization angle maps exhibit similar behavior. 
While complex at low frequency, the polarization angle pat- 
tern becomes more regular in the NGP region at the highest 
frequencies. At 610 MHz, the region of ordered behavior is 
limited to the very high latitude areas. Order expands to 
lower latitudes with increasing frequency, reaching b ~ 40°- 
50° at 1411 MHz. 

It is easy to interpret this in light of the previous dis- 
cussion. The randomization of the polarization angles due to 
Faraday rotation, transferring the power from large to small 
scales, destroys the polarized emission on the largest scales 
at the lowest frequencies. At higher frequencies, the effects of 
Faraday rotation decrease, leading to the re-appearing of the 



large scale structures. This starts in the areas with smaller 
EM. Considering the EM behavior with Galactic latitudes, 
we expect that this would start at the very high latitudes 
and expand to the mid latitudes with increasing frequency, 
producing larger ordered regions showing the intrinsic struc- 
ture. 

This analysis of the iBrouw Spoelstral ()l976f) data 
supports that our patch, given the observing parameters (at 
1.4 GHz and \b\ ^ 40°), is in an intermediate state between 
strong and negligible influence of the Faraday rotation. 



7 IMPLICATIONS FOR CMBP 

The E- and B-mode spectra of Section El are the first mea- 
surements for a high Galactic latitude patch at 1.4 GHz on 
sub-degree scales (i.e. scales on which CMBP has most of 
the power). 

Other spectr a out of the Galac tic plane (b < 30°) have 
been measured bv lHaverkorn et aJ] (|2003h . but at lower fre- 
quency (^ 350 MHz). Their results show slopes with a very 
large distribution, suggesting significant changes of Faraday 
rotation. However, a decrease in slope toward higher Galac- 
tic latitudes seems to exist, in agreement with the interpre- 
tation given here. 

Our data are at a significantly higher frequency and, 
being less affected by Faraday rotation features, are to date 
the most suitable for extrapolation to the frequency range 
of CMB measurements. 

Figure ITI71 shows the ^-mode spectrum extrapolated up 
to the 32 and 90 GHz frequencies of the BaR-SPOrt exper- 
iment, assumed the synchrotron follows a power law 



^synch OC U 



(8) 



with spectral i ndex 7 — —3.1, typ ical of the 1.4-23 GHz 
spectral range ( Bernardi et al. 2004). These spectra include 
the correction for the square of the conversion factor 



2 sinh(x/2) 



X - hu/kTcmh 



(9) 



transforming antenna into thermodynamic temperature of 
CMB (z/ is the frequency and Tcmb = 2.726 K). 

The presence of Faraday rotation effects makes the 
power in our images an upper limit to the intrinsic power 
on CMBP scales. Hence, apart from errors in the frequency 
spectral index, these extrapolations represent an upper limit 
of the contamination on the cosmic signal. 

These results suggest that the synchrotron emission 
should only marginally contaminate the cosmological signal 
at 32 GHz, making this patch a good target for CMBP inves- 
tigations. Even assuming an error of A7 = 0.2 on the power 
law, the contamination would be only a factor 2 stronger, 
which does not change the conclusion. 

The situation at 90 GHz is clearer still: the extrapolated 
spectrum is more than 4 orders of magnitude lower than 
the cosmic signal, leaving the CMBP practically uncontam- 
inated by synchrotron pollution. Including the steepening of 
the synchrotron spe ctral index above 23 GHz observed by 
the WMAP team ( Bennett et alJl2003bl) . the conclusion is 
very robust. 

Such low emission of the Galactic synchtron makes this 
area promising even for the weak B-mode. Its emission has 
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Ip 408 PA 408 




0,0012 ^^^^ 0.60 ^^^^ 180 



Figure 9. Polarized intensity (left) and polarization angle maps (right) formed by interpolating the lBrouw Spoelstral (|l97(f) data. 
The maps correspond to 408, 465, 610, 820, 1411 MHz (from top to bottom). The units are Kelvin and degrees, respectively. The maps, 
in Galactic coordinates centred on the Galactic centre, have been convolved with a 4° FWHM Gaussian filter. 



High Galactic latitude polarized emission at I.4 GHz 11 





10-10 




E 


— m ode 






^ 










10-11 








vAa ; 


cv 












cv 


10-12 


r 


/ 








10-13 






















+ 


10-14 


r 






1 




1 R 












10 












100 






1000 




10-10 




E 


— m ode 












; 




10-11 








vAa ; 


cv 












CV 


10-12 




/ 








10-13 










u 












+ 


10-14 












10-15 








90 GHz , 




10-16 











100 1000 

1 



Figure 10. Fit of the E-mode spectrum measured for our data 
scaled up to 32 (top) and 90 GHz (bottom). The spectrum ex- 
pected for CMBP with cosm olofflcal parameters as from WMAP 
results (f Snergel et ^12003^ is also shown. 
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Figure 11. Fit of the S-mode spectrum measured for our data 
scaled up to 90 GHz. B-mode spectra expected for CMBP with 
tensor-to-scalar perturbations power ratios T/S = 0.1 and T/S = 
0.01 are also shown. The other cosmolo gical parameters are as 
from WMAP results JSpergel et al.ll2003^ . 



a peak near i = 100, whereas our data cover only the i = 
800-2800 range. To compare the emissions, in Figure ITTl we 
extrapolate the synchrotron spectrum down to ^ = 100 with 
the same slope of Table El In fact, since the slope we have 
measured at 1.4 GHz could be steeper than the intrinsic one, 
our estimate is somewhat conservative. 

Figure ITTl indicates that a model with tensor-to-scalar 
perturbation power ratio T/S = 0.1 should be well accessible 
in this area, the synchrotron contribution being a factor 30 
fainter than the CMBP spectrum. This gap makes the result 
quite robust to errors in the extrapolation. 

Even more interesting is that a model with T/S — 0.01 
could be accessible in this part of the sky. At this low level 
it is likely that the leading contaminant will be the thermal 
dust, whose study requires measurements at higher frequen- 
cies (hundreds of GHz). 



8 SUMMARY AND CONCLUSIONS 

We have analysed the observation by Bernardi et al. 
of the polarized emission at 1.4 GHz in the high Galactic 
latitude area {b ^ —40°). This represents the first detection 
of the Galactic synchrotron polarized emission carried out in 
a low emission area at 1.4 GHz on this angular scale, giving 
us good data to estimate the contamination of the CMBP 
signal by the Galactic synchrotron radiation. 

The contamination has been evaluated through the po- 
larized angular power spectra and . These follows a 
power law behavior with spectral indeces (Se = — 1.97±0.08 
and Pe = —1.98 ±0.07. The emission level is about 25 times 
fainter than in Galactic plane regions. 

Extrapolations to the CMB frequency window (30- 
100 GHz) gives encouraging results: the ^-mode of CMBP 
is expected to be safely accessible already at 32 GHz, while 
at 90 GHz the margin is much larger ensuring clean mea- 
surements of CMBP. The low contamination level makes 
this patch a good candidate for the detection of the weaker 
B-mode: the tensorial signal appears accessible for models 
with T/S > 0.01. 

The analysis of the Faraday rotation effects performed 
in this paper reinforces the robustness of our estimates. The 
low BM values measured in the patch (typically 20 rad m~^) 
do not generate significant bandwidth depolarization. Addi- 
tionally, we have carried out a quantitative analysis of the 
effects caused by the randomization of the polarization an- 
gles as a result of Faraday screens. We find an enhancement 
of the power on angular scales smaller than that of random- 
ization, which ensures that the power we measure in the 
observed patch is an upper limit of the intrinsic emission 
- at least on the scales relevant for CMB purposes. More- 
over, the analysis of the Faraday screen effects provides a 
second important result: a steepening of the power law in- 
dex on angular scales smaller than that of randomization. 
The comparison between the slopes we have measured with 
those obtained by other surveys suggests the patch is in an 
intermediate state between significant and negligible Fara- 
day rotation effects. 

A similar conclusion is reach ed from an analysis of the 
data o f iBrouw SDoelstral(ll97(tl . Depolarization effects are 
seen to decrease on the large scale polarized emission as the 
frequency increases from 408 MHz up to 1411 MHz. Start- 
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ing from the NGP at the lowest frequencies, the large scale 
emission appears down to 6 ^ 40° at 1411 MHz. This further 
supports the conclusion that the Galactic latitude of the ob- 
served patch is in an intermediate state between significant 
and negligible Faraday rotation effects at 1.4 GHz. 

It appears that the patch is still affected by Faraday 
effects, but that it is very close to showing the intrinsic 
emission. Therefore, the estimates we have given here are 
conservative upper limits. Observations at higher frequency 
should allow measurements which are not affected by Fara- 
day rotation. 
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